Prediction of m6A and m5C at single-molecule resolution reveals a transcriptome-wide co-occurrence of RNA modifications.

The epitranscriptome embodies many new and largely unexplored functions of RNA. A significant roadblock hindering progress in epitranscriptomics is the identification of more than one modification in individual transcript molecules. We address this with CHEUI (CH3 (methylation) Estimation Using Ionic current). CHEUI predicts N6-methyladenosine (m6A) and 5-methylcytosine (m5C) in individual molecules from the same sample, the stoichiometry at transcript reference sites, and differential methylation between any two conditions. CHEUI processes observed and expected nanopore direct RNA sequencing signals to achieve high single-molecule, transcript-site, and stoichiometry accuracies in multiple tests using synthetic RNA standards and cell line data. CHEUI's capability to identify two modification types in the same sample reveals a co-occurrence of m6A and m5C in individual mRNAs in cell line and tissue transcriptomes. CHEUI provides new avenues to discover and study the function of the epitranscriptome.

Subcellular localization of circular RNAs: Where and why.

Localization of RNAs at specific subcellular locations regulating various local cellular events has gained much attention recently. Like most other classes of RNAs, the function of newly discovered circular RNAs (circRNAs) is predominantly determined by their association with different cellular factors in the cell. CircRNAs function as transcriptional and posttranscriptional regulators of gene expression by interacting with transcription factors, splicing regulators, RNA-binding proteins, and microRNAs or by translating into functional polypeptides. Hence, studying their subcellular localization to assess their function is essential. The discovery of more than a million circRNA and increasing evidence of their involvement in development and diseases require a thorough analysis of their subcellular localization linking to their biological functions. Here, we summarize current knowledge of circRNA localization in cells and extracellular vesicles, factors regulating their subcellular localization, and the implications of circRNA localization on their cellular functions. Given the discovery of many circRNAs in all life forms and their implications in pathophysiology, we discuss the challenges in studying circRNA localization and the opportunities for unlocking the mystery of circRNA functions.

Characterisation of the Arabidopsis thaliana telomerase TERT-TR complex.

Most eukaryotic organisms employ a telomerase complex for the maintenance of chromosome ends. The core of this complex is composed of telomerase reverse transcriptase (TERT) and telomerase RNA (TR) subunits. The TERT reverse transcriptase (RT) domain synthesises telomeric DNA using the TR template sequence. The other TERT domains contribute to this process in different ways. In particular, the TERT RNA-binding domain (TRBD) interacts with specific TR motif(s). Using a yeast 3-hybrid system, we show the critical role of Arabidopsis thaliana (At) TRBD and embryophyta-conserved KRxR motif in the unstructured linker preceding the TRBD domain for binding to the recently identified AtTR subunit. We also show the essential role of the predicted P4 stem and pseudoknot AtTR structures and provide evidence for the binding of AtTRBD to pseudoknot and KRxR motif stabilising interaction with the P4 stem structure. Our results thus provide the first insight into the core part of the plant telomerase complex.

Transfer learning enables identification of multiple types of RNA modifications using nanopore direct RNA sequencing.

Nanopore direct RNA sequencing (DRS) has emerged as a powerful tool for RNA modification identification. However, concurrently detecting multiple types of modifications in a single DRS sample remains a challenge. Here, we develop TandemMod, a transferable deep learning framework capable of detecting multiple types of RNA modifications in single DRS data. To train high-performance TandemMod models, we generate in vitro epitranscriptome datasets from cDNA libraries, containing thousands of transcripts labeled with various types of RNA modifications. We validate the performance of TandemMod on both in vitro transcripts and in vivo human cell lines, confirming its high accuracy for profiling m6A and m5C modification sites. Furthermore, we perform transfer learning for identifying other modifications such as m7G, Ψ, and inosine, significantly reducing training data size and running time without compromising performance. Finally, we apply TandemMod to identify 3 types of RNA modifications in rice grown in different environments, demonstrating its applicability across species and conditions. In summary, we provide a resource with ground-truth labels that can serve as benchmark datasets for nanopore-based modification identification methods, and TandemMod for identifying diverse RNA modifications using a single DRS sample.

CoCoNuTs are a diverse subclass of Type IV restriction systems predicted to target RNA.

A comprehensive census of McrBC systems, among the most common forms of prokaryotic Type IV restriction systems, followed by phylogenetic analysis, reveals their enormous abundance in diverse prokaryotes and a plethora of genomic associations. We focus on a previously uncharacterized branch, which we denote coiled-coil nuclease tandems (CoCoNuTs) for their salient features: the presence of extensive coiled-coil structures and tandem nucleases. The CoCoNuTs alone show extraordinary variety, with three distinct types and multiple subtypes. All CoCoNuTs contain domains predicted to interact with translation system components, such as OB-folds resembling the SmpB protein that binds bacterial transfer-messenger RNA (tmRNA), YTH-like domains that might recognize methylated tmRNA, tRNA, or rRNA, and RNA-binding Hsp70 chaperone homologs, along with RNases, such as HEPN domains, all suggesting that the CoCoNuTs target RNA. Many CoCoNuTs might additionally target DNA, via McrC nuclease homologs. Additional restriction systems, such as Type I RM, BREX, and Druantia Type III, are frequently encoded in the same predicted superoperons. In many of these superoperons, CoCoNuTs are likely regulated by cyclic nucleotides, possibly, RNA fragments with cyclic termini, that bind associated CARF (CRISPR-Associated Rossmann Fold) domains. We hypothesize that the CoCoNuTs, together with the ancillary restriction factors, employ an echeloned defense strategy analogous to that of Type III CRISPR-Cas systems, in which an immune response eliminating virus DNA and/or RNA is launched first, but then, if it fails, an abortive infection response leading to PCD/dormancy via host RNA cleavage takes over.

All organisms, from animals to bacteria, are subject to genetic parasites, such as viruses and transposons. Genetic parasites are pieces of nucleic acids (DNA or RNA) that can use a cell's machinery to copy themselves at the expense of their hosts. This often leads to the host's demise, so organisms evolved many types of defense mechanisms. One of the most ancient and common forms of defense against viruses and transposons is the targeted restriction of nucleic acids, that is, deployment of host enzymes that can destroy or restrict nucleic acids containing specific sequence motifs or modifications. In bacteria, many of the restriction enzymes targeting parasitic genetic elements are formed by fusions of proteins from the so-called McrBC systems with a protein domain called EVE. EVE and other functionally similar domains are a part of proteins that recognize and bind modified bases in nucleic acids. Enzymes can use the ability of these specificity domains to bind modified bases to detect non-host nucleic acids. Bell et al. conducted a comprehensive computational search for McrBC systems and discovered a large and highly diverse branch of this family with unusual characteristic structural and functional domains. These features include regions that form long alpha-helices (coils) that coil with other alpha-helices (known as coiled-coils), as well as several distinct enzymatic domains that break down nucleic acids (known as nucleases). They call these systems CoCoNuTs (coiled-coiled nuclease tandems). All CoCoNuTs contain domains, including EVE-like ones, which are predicted to interact with components of the RNA-based systems responsible for producing proteins in the cell (translation), suggesting that the CoCoNuTs have an important impact on protein abundance and RNA metabolism. Bell et al.'s findings will be of interest to scientists working on prokaryotic immunity and virulence. Furthermore, similarities between CoCoNuTs and components of eukaryotic RNA-degrading systems suggest evolutionary connections between this diverse family of bacterial predicted RNA restriction systems and RNA regulatory pathways of eukaryotes. Further deciphering the mechanisms of CoCoNuTs could shed light on how certain pathways of RNA metabolism and regulation evolved, and how they may contribute to advances in biotechnology.

Harnessing noncanonical crRNA for highly efficient genome editing.

The CRISPR-Cas12a system is more advantageous than the widely used CRISPR-Cas9 system in terms of specificity and multiplexibility. However, its on-target editing efficiency is typically much lower than that of the CRISPR-Cas9 system. Here we improved its on-target editing efficiency by simply incorporating 2-aminoadenine (base Z, which alters canonical Watson-Crick base pairing) into the crRNA to increase the binding affinity between crRNA and its complementary DNA target. The resulting CRISPR-Cas12a (named zCRISPR-Cas12a thereafter) shows an on-target editing efficiency comparable to that of the CRISPR-Cas9 system but with much lower off-target effects than the CRISPR-Cas9 system in mammalian cells. In addition, zCRISPR-Cas12a can be used for precise gene knock-in and highly efficient multiplex genome editing. Overall, the zCRISPR-Cas12a system is superior to the CRISPR-Cas9 system, and our simple crRNA engineering strategy may be extended to other CRISPR-Cas family members as well as their derivatives.

A comprehensive analysis of spermatozoal RNA elements in idiopathic infertile males undergoing fertility treatment.

Current approaches to diagnosing male infertility inadequately assess the complexity of the male gamete. Beyond the paternal haploid genome, spermatozoa also deliver coding and non-coding RNAs to the oocyte. While sperm-borne RNAs have demonstrated potential involvement in embryo development, the underlying mechanisms remain unclear. In this study, 47 sperm samples from normozoospermic males undergoing fertility treatment using donor oocytes were sequenced and analyzed to evaluate associations between sperm RNA elements (exon-sized sequences) and blastocyst progression. A total of 366 RNA elements (REs) were significantly associated with blastocyst rate (padj < 0.05), some of which were linked to genes related to critical developmental processes, including mitotic spindle formation and both ectoderm and mesoderm specification. Of note, 27 RE-associated RNAs are predicted targets of our previously reported list of developmentally significant miRNAs. Inverse RE-miRNA expression patterns were consistent with miRNA-mediated down-regulation. This study provides a comprehensive set of REs which differ by the patient's ability to produce blastocysts. This knowledge can be leveraged to improve clinical screening of male infertility and ultimately reduce time to pregnancy.

TransAC4C-a novel interpretable architecture for multi-species identification of N4-acetylcytidine sites in RNA with single-base resolution.

N4-acetylcytidine (ac4C) is a modification found in ribonucleic acid (RNA) related to diseases. Expensive and labor-intensive methods hindered the exploration of ac4C mechanisms and the development of specific anti-ac4C drugs. Therefore, an advanced prediction model for ac4C in RNA is urgently needed. Despite the construction of various prediction models, several limitations exist: (1) insufficient resolution at base level for ac4C sites; (2) lack of information on species other than Homo sapiens; (3) lack of information on RNA other than mRNA; and (4) lack of interpretation for each prediction. In light of these limitations, we have reconstructed the previous benchmark dataset and introduced a new dataset including balanced RNA sequences from multiple species and RNA types, while also providing base-level resolution for ac4C sites. Additionally, we have proposed a novel transformer-based architecture and pipeline for predicting ac4C sites, allowing for highly accurate predictions, visually interpretable results and no restrictions on the length of input RNA sequences. Statistically, our work has improved the accuracy of predicting specific ac4C sites in multiple species from less than 40% to around 85%, achieving a high AUC > 0.9. These results significantly surpass the performance of all existing models.

A comparison of five methods to maximize RNA and DNA isolation yield from adipose tissue.

Adipose tissue in the human body occurs in various forms with different functions. It is an energy store, a complex endocrine organ, and a source of cells used in medicine. Many molecular analyses require the isolation of nucleic acids, which can cause some difficulties connected with the large amount of lipids in adipocytes. Ribonucleic acid isolation is particularly challenging due to its low stability and easy degradation by ribonucleases. The study aimed to compare and evaluate five RNA and DNA isolation methods from adipose tissue. The tested material was subcutaneous porcine adipose tissue subjected to different homogenization methods and RNA or DNA purification. A mortar and liquid nitrogen or ceramic beads were used for homogenization. The organic extraction (TriPure Reagent), spin columns with silica-membrane (RNeasy Mini Kit or High Pure PCR Template Preparation Kit), and the automatic MagNA Pure system were used for the purification. Five combinations were compared for RNA and DNA isolation. Obtained samples were evaluated for quantity and quality. The methods were compared in terms of yield (according to tissue mass), purity (A260/280 and A260/230), and nucleic acid degradation (RNA Integrity Number, RIN; DNA Integrity Number, DIN). The results were analyzed statistically. The average RNA yield was highest in method I, which used homogenization with ceramic beads and organic extraction. Low RNA concentration didn't allow us to measure degradation for all samples in method III (homogenization with ceramic beads and spin-column purification). The highest RNA quality was achieved with method IV using homogenization in liquid nitrogen and spin column purification, which makes it the most effective for RNA isolation from adipose tissue. Required values of DNA yield, purity, and integrity were achieved only with spin column-based methods (III and IV). The most effective method for DNA isolation from adipose tissue is method III, using spin-columns without additional homogenization.

Noncanonical functions of telomerase and telomeres in viruses-associated cancer.

The cause of cancer is attributed to the uncontrolled growth and proliferation of cells resulting from genetic changes and alterations in cell behavior, a phenomenon known as epigenetics. Telomeres, protective caps on the ends of chromosomes, regulate both cellular aging and cancer formation. In most cancers, telomerase is upregulated, with the telomerase reverse transcriptase (TERT) enzyme and telomerase RNA component (TERC) RNA element contributing to the maintenance of telomere length. Additionally, it is noteworthy that two viruses, human papillomavirus (HPV) and Epstein-Barr virus (EBV), utilize telomerase for their replication or persistence in infected cells. Also, TERT and TERC may play major roles in cancer not related to telomere biology. They are involved in the regulation of gene expression, signal transduction pathways, cellular metabolism, or even immune response modulation. Furthermore, the crosstalk between TERT, TERC, RNA-binding proteins, and microRNAs contributes to a greater extent to cancer biology. To understand the multifaceted roles played by TERT and TERC in cancer and viral life cycles, and then to develop effective therapeutic strategies against these diseases, are fundamental for this goal. By investigating deeply, the complicated mechanisms and relationships between TERT and TERC, scientists will open the doors to new therapies. In its analysis, the review emphasizes the significance of gaining insight into the multifaceted roles that TERT and TERC play in cancer pathogenesis, as well as their involvement in the viral life cycle for designing effective anticancer therapy approaches.

Host RNA Expression Signatures in Young Infants with Urinary Tract Infection: A Prospective Study.

Early diagnosis of infections in young infants remains a clinical challenge. Young infants are particularly vulnerable to infection, and it is often difficult to clinically distinguish between bacterial and viral infections. Urinary tract infection (UTI) is the most common bacterial infection in young infants, and the incidence of associated bacteremia has decreased in the recent decades. Host RNA expression signatures have shown great promise for distinguishing bacterial from viral infections in young infants. This prospective study included 121 young infants admitted to four pediatric emergency care departments in the capital region of Denmark due to symptoms of infection. We collected whole blood samples and performed differential gene expression analysis. Further, we tested the classification performance of a two-gene host RNA expression signature approaching clinical implementation. Several genes were differentially expressed between young infants with UTI without bacteremia and viral infection. However, limited immunological response was detected in UTI without bacteremia compared to a more pronounced response in viral infection. The performance of the two-gene signature was limited, especially in cases of UTI without bloodstream involvement. Our results indicate a need for further investigation and consideration of UTI in young infants before implementing host RNA expression signatures in clinical practice.

Unearthing a novel function of SRSF1 in binding and unfolding of RNA G-quadruplexes.

SRSF1 governs splicing of over 1500 mRNA transcripts. SRSF1 contains two RNA-recognition motifs (RRMs) and a C-terminal Arg/Ser-rich region (RS). It has been thought that SRSF1 RRMs exclusively recognize single-stranded exonic splicing enhancers, while RS lacks RNA-binding specificity. With our success in solving the insolubility problem of SRSF1, we can explore the unknown RNA-binding landscape of SRSF1. We find that SRSF1 RS prefers purine over pyrimidine. Moreover, SRSF1 binds to the G-quadruplex (GQ) from the ARPC2 mRNA, with both RRMs and RS being crucial. Our binding assays show that the traditional RNA-binding sites on the RRM tandem and the Arg in RS are responsible for GQ binding. Interestingly, our FRET and circular dichroism data reveal that SRSF1 unfolds the ARPC2 GQ, with RS leading unfolding and RRMs aiding. Our saturation transfer difference NMR results discover that Arg residues in SRSF1 RS interact with the guanine base but not other nucleobases, underscoring the uniqueness of the Arg/guanine interaction. Our luciferase assays confirm that SRSF1 can alleviate the inhibitory effect of GQ on gene expression in the cell. Given the prevalence of RNA GQ and SR proteins, our findings unveil unexplored SR protein functions with broad implications in RNA splicing and translation.

Discovery of a polymorphic gene fusion via bottom-up chimeric RNA prediction.

Gene fusions and their chimeric products are commonly linked with cancer. However, recent studies have found chimeric transcripts in non-cancer tissues and cell lines. Large-scale efforts to annotate structural variations have identified gene fusions capable of generating chimeric transcripts even in normal tissues. In this study, we present a bottom-up approach targeting population-specific chimeric RNAs, identifying 58 such instances in the GTEx cohort, including notable cases such as SUZ12P1-CRLF3, TFG-ADGRG7 and TRPM4-PPFIA3, which possess distinct patterns across different ancestry groups. We provide direct evidence for an additional 29 polymorphic chimeric RNAs with associated structural variants, revealing 13 novel rare structural variants. Additionally, we utilize the All of Us dataset and a large cohort of clinical samples to characterize the association of the SUZ12P1-CRLF3-causing variant with patient phenotypes. Our study showcases SUZ12P1-CRLF3 as a representative example, illustrating the identification of elusive structural variants by focusing on those producing population-specific fusion transcripts.

MazF-rolling circle amplification combined MALDI-TOF MS for site-specific detection of N6-methyladenosine RNA.

N6-methyladenosine (m6A) is one of the most abundant chemical modifications in RNA and has vital significance in cellular processes and tumor development. However, the accurate analysis of site-specific m6A modification remains a challenge. In this work, a MazF endoribonuclease activated rolling circle amplification (MazF-RCA) combined MALDI-TOF MS assay is developed for the detection of site-specific m6A-RNA. MazF endoribonuclease can specifically cleave the ACA motif, leaving methylated (m6A)CA motif intact. The intact methylated RNA can then be amplified through rolling circle amplification, and the generated reporter oligonucleotides are detected by MALDI-TOF MS. The assay exhibits good quantification ability, presenting a wide linear range (100 fM to 10 nM) with the limit-of-detection lower than 100 fM. Additionally, the assay can accurately detect methylated RNA in the presence of large amount of non-methylated RNA with a relative abundance of methylated RNA down to 0.5%. The developed assay was further applied to detect m6A-RNA spiked in MCF-7 cell RNA extracts, with the recovery rates in the range of 90.64-106.93%. The present assay provides a novel platform for the analysis of site-specific m6A-RNA at high specificity and sensitivity, which can promote the study of RNA methylation in clinical and biomedical research.

Detection of newly synthesized RNA reveals transcriptional reprogramming during ZGA and a role of Obox3 in totipotency acquisition.

Zygotic genome activation (ZGA) after fertilization enables the maternal-to-zygotic transition. However, the global view of ZGA, particularly at initiation, is incompletely understood. Here, we develop a method to capture and sequence newly synthesized RNA in early mouse embryos, providing a view of transcriptional reprogramming during ZGA. Our data demonstrate that major ZGA gene activation begins earlier than previously thought. Furthermore, we identify a set of genes activated during minor ZGA, the promoters of which show enrichment of the Obox factor motif, and find that Obox3 or Obox5 overexpression in mouse embryonic stem cells activates ZGA genes. Notably, the expression of Obox factors is severely impaired in somatic cell nuclear transfer (SCNT) embryos, and restoration of Obox3 expression corrects the ZGA profile and greatly improves SCNT embryo development. Hence, our study reveals dynamic transcriptional reprogramming during ZGA and underscores the crucial role of Obox3 in facilitating totipotency acquisition.

Clinical mutations in the TERT and TERC genes coding for telomerase components induced oxidative stress, DNA damage at telomeres and cell apoptosis besides decreased telomerase activity.

Telomeres are nucleoprotein structures at the end of chromosomes that maintain their integrity. Mutations in genes coding for proteins involved in telomere protection and elongation produce diseases such as dyskeratosis congenita or idiopathic pulmonary fibrosis known as telomeropathies. These diseases are characterized by premature telomere shortening, increased DNA damage and oxidative stress. Genetic diagnosis of telomeropathy patients has identified mutations in the genes TERT and TERC coding for telomerase components but the functional consequences of many of these mutations still have to be experimentally demonstrated. The activity of twelve TERT and five TERC mutants, five of them identified in Spanish patients, has been analyzed. TERT and TERC mutants were expressed in VA-13 human cells that express low telomerase levels and the activity induced was analyzed. The production of reactive oxygen species, DNA oxidation and TRF2 association at telomeres, DNA damage response and cell apoptosis were determined. Most mutations presented decreased telomerase activity, as compared to wild-type TERT and TERC. In addition, the expression of several TERT and TERC mutants induced oxidative stress, DNA oxidation, DNA damage, decreased recruitment of the shelterin component TRF2 to telomeres and increased apoptosis. These observations might indicate that the increase in DNA damage and oxidative stress observed in cells from telomeropathy patients is dependent on their TERT or TERC mutations. Therefore, analysis of the effect of TERT and TERC mutations of unknown function on DNA damage and oxidative stress could be of great utility to determine the possible pathogenicity of these variants.

Detecting m6A at single-molecular resolution via direct RNA sequencing and realistic training data.

Direct RNA sequencing offers the possibility to simultaneously identify canonical bases and epi-transcriptomic modifications in each single RNA molecule. Thus far, the development of computational methods has been hampered by the lack of biologically realistic training data that carries modification labels at molecular resolution. Here, we report on the synthesis of such samples and the development of a bespoke algorithm, mAFiA (m6A Finding Algorithm), that accurately detects single m6A nucleotides in both synthetic RNAs and natural mRNA on single read level. Our approach uncovers distinct modification patterns in single molecules that would appear identical at the ensemble level. Compared to existing methods, mAFiA also demonstrates improved accuracy in measuring site-level m6A stoichiometry in biological samples.

The dynamic landscape of enhancer-derived RNA during mouse early embryo development.

Enhancer-derived RNAs (eRNAs) play critical roles in diverse biological processes by facilitating their target gene expression. However, the abundance and function of eRNAs in early embryos are not clear. Here, we present a comprehensive eRNA atlas by systematically integrating publicly available datasets of mouse early embryos. We characterize the transcriptional and regulatory network of eRNAs and show that different embryo developmental stages have distinct eRNA expression and regulatory profiles. Paternal eRNAs are activated asymmetrically during zygotic genome activation (ZGA). Moreover, we identify an eRNA, MZGAe1, which plays an important function in regulating mouse ZGA and early embryo development. MZGAe1 knockdown leads to a developmental block from 2-cell embryo to blastocyst. We create an online data portal, M2ED2, to query and visualize eRNA expression and regulation. Our study thus provides a systematic landscape of eRNA and reveals the important role of eRNAs in regulating mouse early embryo development.

BEERS2: RNA-Seq simulation through high fidelity in silico modeling.

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully length messenger RNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in polymerase chain reaction (PCR) amplification, barcode read errors and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

Circular RNAs: Biogenesis, Functions, and Role in Myocardial Hypertrophy.

Circular RNAs (circRNAs) are a large class of endogenous single-stranded covalently closed RNA molecules. High-throughput RNA sequencing and bioinformatic algorithms have identified thousands of eukaryotic circRNAs characterized by high stability and tissue-specific expression pattern. Recent studies have shown that circRNAs play an important role in the regulation of physiological processes in the norm and in various diseases, including cardiovascular disorders. The review presents current concepts of circRNA biogenesis, structural features, and biological functions, describes the methods of circRNA analysis, and summarizes the results of studies on the role of circRNAs in the pathogenesis of hypertrophic cardiomyopathy, the most common inherited heart disease.